Technical  Report  W-94-4 
November  1 994 


US  Army  Corps 
of  Engineers 


Waterways  Experiment 


Station 


Water  Quality  Research  Program 

RECOVERY,  A  Mathematical  Model 
to  Predict  the  Temporal  Response 
of  Surface  Water  to  Contaminated  Sediments 

by  Jean  M.  Boyer,  Steven  C.  Chapra,  University  of  Colorado 
Carlos  E.  Ruiz,  Mark  S.  Dortch,  WES 


Approved  For  Public  Release;  Distribution  Is  Unlimited 


Prepared  for  Headquarters,  U.S.  Army  Corps  of  Engineers 
U.S.  Army  Engineer  Division,  New  England 
and  U.S.  Army  Engineer  District,  New  York 


The  contents  of  this  report  are  not  to  be  used  for  advertising, 
publication,  or  promotional  purposes.  Citation  of  tradenames 
does  not  constitute  an  official  endorsement  or  approval  of  the  use 
of  such  commercial  products. 


PRINTED  ON  RECYCLED  PAPER 


Water  Quality 
Research  Program 


Technical  Report  W-94-4 
November  1 994 


RECOVERY,  A  Mathematical  Model 
to  Predict  the  Temporal  Response 
of  Surface  Water  to  Contaminated  Sediments 

by  Jean  M.  Boyer,  Steven  C.  Chapra 

University  of  Colorado 
Center  for  Advanced  Decision  Support 
for  Water  and  Environmental  Systems 
Boulder,  CO  80309-0428 

Carlos  E.  Ruiz,  Mark  S.  Dortch 

U.S.  Army  Corps  of  Engineers 
Waterways  Experiment  Station 
3909  Halls  Ferry  Road 
Vicksburg,  MS  39180-6199 


twto  ott.at.itt  inspected  q 

Final  report  i'" 

Approved  for  public  release;  distribution  is  unlimited 


Prepared  for  U.S.  Army  Corps  of  Engineers 
Washington,  DC  20314-1000 
U.S.  Army  Engineer  Division,  New  England 
Waltham,  MA  02254-9149 
U.S.  Army  Engineer  District,  New  York 
New  York,  NY  10278-0090 


and 


Waterways  Experiment  Station  Cataloging-in-Publication  Data 

RECOVERY :  a  mathematical  model  to  predict  the  temporal  response  of 
surface  water  to  contaminated  sediments  /  by  Jean  M.  Boyer ...  [et  al.] ; 
prepared  for  U.S.  Army  Corps  of  Engineers,  U.S.  Army  Engineer  Divi¬ 
sion,  New  England  and  U.S.  Army  Engineer  District,  New  York. 

61  p. :  ill. ;  28  cm.  —  (Technical  report ;  W-94-4) 

Includes  bibliographic  references. 

1 .  Water  quality  —  Measurement  —  Data  processing.  2.  RECOVERY 
(Computer  program)  3.  Marine  sediments  —  Environmental  aspects. 

4.  Sediments  (Geology)  —  Environmental  aspects.  I.  Boyer,  Jean  M. 

II.  United  States.  Army.  Corps  of  Engineers.  III.  United  States.  Army. 
Corps  of  Engineers.  New  England  Division.  IV.  United  States.  Army. 
Corps  of  Engineers.  New  York  District.  V.  U.S.  Army  Engineer  Waterways 
Experiment  Station.  VI.  Water  Quality  Research  Program.  VII.  Title:  A 
mathematical  model  to  predict  the  temporal  response  of  surface  water  to 
contaminated  sediments.  VIII.  Series:  Technical  report  (U.S.  Army  Engi¬ 
neer  Waterways  Experiment  Station) ;  W-94-4. 

TA7  W34  no.W-94-4 


Preface . v 

1 —  Introduction . 1 

2 —  Model  Description  . 2 

Contaminant  Mass  Balances . 2 

Solids  Budget . 7 

Default  Parameters  . 8 

Numerical  Methods .  9 

3 —  Model  Confirmation  . . 10 

Study  Area .  10 

Overview .  11 

Rate  Coefficients .  13 

Simulation  Results .  14 

4 —  User’s  Guide  . 20 

Changing  Chemical  Database . 22 

Model  Source . 22 

Program  Features  . 23 

System  Specific  Issues . 26 

5 —  Conclusions  and  Recommendations  . 28 

References  . 30 

Appendix  A:  Nomenclature . A1 

Appendix  B:  Input  Screens  for  DDE  and  Lindane  Runs  . B1 

Appendix  C:  Screens  . Cl 

SF  298 


List  of  Figures 


Figure  1.  Physical  configuration  of  RECOVERY 
Figure  2.  Schematic  of  RECOVERY . 


2 

3 


Figure  3.  Lindane  in  water  column  . 15 

Figure  4.  Lindane  in  sediment . 15 

Figure  5.  Volatilization  rate  of  lindane  from  water 

computed  by  model  . 16 

Figure  6.  Flux  of  lindane  from  sediments  into  water 

computed  by  model  . 16 

Figure  7.  DDE  in  water  column . 17 

Figure  8.  DDE  in  sediment . 17 

Figure  9.  Volatilization  rate  of  DDE  from  water  computed 

by  model  . 18 

Figure  10.  Flux  of  DDE  from  sediments  into  water 

computed  by  model  . 18 

Figure  11.  Final  simulation  of  DDE  in  water  column . 19 

Figure  12.  Final  simulation  of  DDE  in  sediment . 19 


List  of  Tables 


Table  1.  List  of  Assumptions  Used  in  Development  of  RECOVERY  ...  3 


Table  2.  Concentrations  of  DDE  and  Lindane  in  Water  Column . 11 

Table  3.  DDE  and  Lindane  Concentration  in  Bottom  Sediments . 12 

Table  4.  Input  Data  for  RECOVERY  Model  . 24 

Table  5.  Compound  Properties . 25 


Preface 


The  work  reported  herein  was  sponsored  by  the  U.S.  Army  Engineer 
Division,  New  England,  U.S.  Army  Engineer  District,  New  York,  and  by  the 
Headquarters,  U.S.  Army  Corps  of  Engineers  (HQUSACE),  as  part  of  the 
Water  Quality  Research  Program  (WQRP),  Work  Unit  32808.  The  WQRP  is 
assigned  to  the  U.S.  Army  Engineer  Waterways  Experiment  Station  (WES) 
under  the  purview  of  the  Environmental  Laboratory  (EL).  Funding  was 
provided  under  Department  of  the  Army  Appropriation  No.  96X3121,  General 
Investigation.  The  WQRP  is  managed  under  the  Environmental  Resources 
Research  and  Assistance  Programs  (ERRAP),  Mr.  J.  L.  Decell,  Manager. 

Mr.  Robert  C.  Gunkel  was  Assistant  Manger,  ERRAP,  for  the  WQRP.  Tech¬ 
nical  Monitors  during  this  study  were  Mr.  Frederick  B.  Juhle,  Mr.  Rixie 
Hardy,  and  Dr.  John  Bushman,  HQUSACE. 

This  report  provides  documentation  of  the  RECOVERY  model 
(Version  1.0),  including  mathematical  formulations,  computational  algorithms, 
a  user’s  guide,  and  a  confirmation  application.  The  Principal  Investigator  for 
this  study  was  Dr.  Carlos  E.  Ruiz,  Water  Quality  and  Contaminant  Modeling 
Branch  (WQCMB),  Environmental  Processes  and  Effects  Division  (EPED), 
EL.  Model  development  was  conducted  under  contract  to  the  Center  for 
Advanced  Decision  Support  for  Water  and  Environmental  Systems 
(CADSWES),  University  of  Colorado,  Boulder,  CO.  Dr.  Steven  C.  Chapra 
and  Ms.  Jean  M.  Boyer  of  CADSWES  performed  the  original  model  develop¬ 
ment.  Drs.  Ruiz  and  Marie  S.  Dortch  of  WQCMB  provided  guidance  during 
development. 

This  report  was  prepared  by  Ms.  Boyer  and  Drs.  Chapra,  Ruiz,  and  Dortch. 
This  study  was  conducted  under  the  direct  supervision  of  Dr.  Dortch,  Chief, 
WQCMB,  and  under  the  general  supervision  of  Mr.  Donald  L.  Robey,  Chief, 
EPED,  and  Dr.  John  W.  Keeley,  Director,  EL.  This  report  was  reviewed  by 
Drs.  Douglas  Gunnison,  Environmental  Processes  and  Effects  Branch,  EPED, 
and  Carl  F.  Cerco,  WQCMB. 

At  the  time  of  publication  of  this  report,  Director  of  WES  was 
Dr.  Robert  W.  Whalin.  Commander  was  COL  Bruce  K.  Howard,  EN. 


v 


This  report  should  be  cited  as  follows: 

Boyer,  J.  M.,  Chapra,  S.  C.,  Ruiz,  C.  E.,  and  Dortch,  M.  S. 
(1994).  “RECOVERY,  a  mathematical  model  to  predict  the 
temporal  response  of  surface  water  to  contaminated  sediments,” 
Technical  Report  W-94-4,  U.S.  Army  Engineer  Waterways 
Experiment  Station,  Vicksburg,  MS. 


The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication , 
or  promotional  purposes.  Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of  such  commercial  products. 


1  Introduction 


Toxic  substances  are  known  to  associate  strongly  with  fine-grain  particulate 
matter.  Consequently,  it  is  not  surprising  that  toxicants  tend  to  accumulate  in 
the  bottom  sediments  of  natural  waters.  An  important  question  in  environmen¬ 
tal  impact  assessment  is  to  what  extent  do  these  contaminated  sediments  influ¬ 
ence  overlying  waters. 

The  report  herein  describes  a  modeling  framework  that  is  designed  to  assess 
the  impact  of  contaminated  bottom  sediments  on  surface  waters.  The  analysis 
is  limited  to  cases  where  the  overlying  water  is  well  mixed.  In  addition,  the 
contaminant  is  assumed  to  be  organic  and  to  follow  linear  equilibrium  sorption 
and  first-order  decay  kinetics. 

The  framework  is  implemented  via  a  personal  computer  software 
package — RECOVERY.  It  is  designed  for  interactive  implementation  and  is 
self-documented.  RECOVERY  allows  the  user  to  rapidly  generate  and  analyze 
recovery  scenarios  for  contaminated  systems.  The  software  includes  graphical 
displays  and  can  be  run  on  IBM-PC  compatible  microcomputers. 

This  report  provides  supplemental  documentation,  including  a  mathematical 
description  of  the  model,  its  computational  algorithm,  and  a  confirmation 
implementation.  It  also  includes  a  user’s  guide. 
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2  Model  Description 


The  present  model  is  an  application  and  extension  of  frameworks  developed 
previously  (Chapra  1982,  1986;  Chapra  and  Reckhow  1983).  As  in  Figures  1 
and  2,  the  system  is  idealized  as  a  well-mixed  surface  water  layer  underlain  by 
a  vertically  stratified  sediment  column.  The  sediment  is  well  mixed  horizon¬ 
tally,  but  segmented  vertically  into  a  well-mixed  surface  layer  and  a  deep  sedi¬ 
ment.  The  latter,  in  turn,  is  segmented  into  contaminated  and  uncontaminated 
regions.  The  specification  of  a  mixed  surface  layer  is  included  because  an 
unconsolidated  layer  is  often  observed  at  the  surface  of  sediments  because  of  a 
number  of  processes,  including  bioturbation  and  mechanical  mixing. 

A  list  of  assumptions 
used  in  the  development  of 
this  model  can  be  found  in 
Table  1.  Also,  nomenclature 
used  in  this  report  is  listed 
in  Appendix  A. 
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Figure  1 .  Physical  configuration  of 
RECOVERY 

Contaminant  Mass  Balances 

A  mass  balance  for  the  contaminant  in  the  water  can  be  written  as 


dc 

V _  =  Qc  -  Qc  -  k  V  c  -  kV  c  -  v  A  F  c  +  v  A  c 

w  ^  i  w  www  v  w  w  s  w  pw  w  r  m  n 


+  VA(Vm  '  F*»C») 


(1) 
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Figure  2.  Schematic  ot  RECOVERY 


Table  1 

List  of  Assumptions  Used  in  Development  of  RECOVERY 

Contaminant  is  of  organic  nature. 

Water  body  is  well  mixed. 

Mixed  layer  is  well  mixed. 

In  the  deep  sediments,  contaminant  concentration  varies  in  vertical  direction  only. 

Initial  concentration  of  compound  in  contaminated  region  is  uniform  throughout  region. 

Initial  concentration  of  contaminant  in  region  below  contaminated  region  is  zero. 

Sediments  are  only  source  of  contaminant  to  water  body. 

Contaminant  follows  a  linear  equilibrium  sorption  mechanism. 

Contaminant  decays  according  to  first  order  kinetics. 

There  is  no  compaction  in  sediments. 

System  is  isothermal. 

Water  body  is  at  25  °C. 

Only  single  component  mass  transfer  is  described.  Movement  of  contaminant  is  independent 
of  presence  of  other  contaminants. 

Flow  through  is  constant. 
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where 


Vw  =  volume  of  water  body,  m3 

cw  and  cm  =  concentrations  of  toxicant  in  water  and  mixed  sediments, 
respectively,  pg/m3 

t  -  time,  years 

Q  =  flushing  flow  rate,  m3/year 

c,  =  inflow  concentration,  pg/m3,  which  reflects  both  direct  and 
tributary  loadings  (in  the  present  version  of  the  model,  this 
concentration  is  assumed  to  be  zero) 

kw  =  decay  rate  constant  of  contaminant  in  water,  /year 

kv  =  volatilization  rate  of  contaminant,  /year 

vs  =  settling  velocity  of  particulate  matter,  m/year 

Aw  and  Am  =  surface  areas  of  water  and  mixed  sediment,  respectively,  m2 

=  fraction  of  contaminant  in  particulate  form  in  water 

vr  =  resuspension  velocity  of  sediments,  m/year 

vd  =  diffusion  mass  transfer  coefficient  at  sediment-water 
interface,  m/year 

Fjp  =  ratio  of  contaminant  concentration  in  sediment  pore  water 
to  contaminant  concentration  in  total  sediment 

=  fraction  of  contaminant  in  dissolved  form  in  water 


The  initial  condition  for  Equation  1  at  t  =  0  is  cw  =  cwo. 

A  mass  balance  for  the  mixed-sediment  layer  can  be  written  as 

dc 

VnT-T-  =  +  ^  AJ  CyJ  ~ 

dt  (2) 

+  +  v/JF+c,( 0)  - 
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where 


Vm  -  volume  of  mixed  layer,  m3 

km  =  decay  rate  constant  of  contaminant  in  mixed  layer,  /year 
vh  =  burial  velocity,  m/year 

cJO)  =  contaminant  concentration  at  top  of  deep  contaminated 
layer,  pg/m3 

The  initial  condition  for  this  equation  at  t  =  0  is  cm  =  c^. 

Both  the  deep  contaminated  and  clean  sediments  can  be  modeled  with  one¬ 
dimensional  advection-diffusion-decay  equations  of  the  form 


~3F 


-  kc. 


(3) 


where 

cs  =  contaminant  concentration  in  deep  sediments,  pg/m3 
<J)  =  the  sediment  porosity 

Ds  =  diffusion  rate  in  sediment  pore  water,  m2/year 
i  -  depth  into  sediment,  m,  where  z  =  0  at  top  of  deep  sediments 
ks  -  decay  rate  constant  of  contaminant  in  deep  sediments,  /year 

The  initial  and  boundary  conditions  for  Equation  3  are  as  follows: 


IC: 

at 

t  =  0,  cs  =  cso  ( zm<z<L ) 

IC: 

at 

t  =  0,  cs  =  0  (L<z<«>) 

BC  1: 

at 

^  =  ^m>  3  ~  J ms 

BC  2: 

at 

dc 

2  =  00  ’  ^  =  0 
az 
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where 


L  =  distance  from  top  of  mixed  layer  to  bottom  of  contaminated  layer 

J  =  mass  flux  of  contaminant,  g/(m2-year) 

Jm  =  mass  flux  of  contaminant  from  mixed  layer  to  sediment  layer, 
g/(m2-year) 

The  F's  in  Equations  1-3  are  dimensionless  ratios  that  reflect  how  the  con¬ 
taminant  partitions  between  solid  matter  and  water  assuming  a  linear  sorption 
mechanism.  Denoting  the  contaminant  as  A,  they  are  defined  as 


Mass  of  A  ( Particulate  Form) 
Total  Volume 


Mass  of  A  ( Dissolved  Form)  +  Mass  of  A  (Particulate  Form) 
Total  Volume 


(4) 


JWattr  body 


F '  -  1  -  F  =  1 

1  +  S 


Mass  of  A  [Dissolved  Form ) 


Total  Volume 


Mass  of  A  [Dissolved  Form)  +  Mass  of  A  [Particulate  Form) 


Total  Volume 


(5) 


JWater  Body 


*  4>  +  A*  (l  -  <>)p 


/ 

Mass  of  A  {Dissolved  Form) 

Volume  of  Pore  Water 

Mass  of  A  ( Dissolved  Form)  +  Mass  of  A  ( Particulate  Form) 
Total  Volume 

\ 


jSediment 


(6) 


where 

**  and  =  contaminant  partition  coefficient  in  water  and  sediments, 

respectively,  m3/g 

p  =  density  of  sediment  solids,  g/m3 
sw  =  suspended  solids  concentration  in  water,  g/m3 
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Note  that  the  model  allows  different  values  of  this  parameter  to  be  employed 
in  the  mixed  layer  and  the  vertically  segmented  sediments  (F<ipm  and  F^j. 

The  decay  rate  constants  km,  and  ks  represent  all  mechanisms  for  decay 
except  volatilization,  which  is  accounted  for  separately.  These  mechanisms 
include  photolysis,  hydrolysis,  and  bacterial  degradation. 

The  mass  transfer  coefficient  for  diffusive  sediment-water  exchange  is 
related  to  more  fundamental  parameters  by 


where 

(j)  =  porosity 

D3  =  diffusion  coefficient  in  sediment  pore  water,  m2/year 

z'  =  characteristic  length  over  which  gradient  exists  at  sediment- 
water  interface,  m. 

Note  that  a  value  of  1  cm  is  assumed  for  z'  after  Thomann  and  Mueller 
(1987).  Also,  Ds  is  related  to  molecular  diffusivity  Dm  by  the  relation  (Berner 
1980;  Manheim  and  Waterman  1974) 

D.  =  DJ2 


Solids  Budget 

The  velocity  terms — v„  vr,  and  vt — in  Equations  1  and  2  are  computed 
according  to  a  steady-state  mass  balance  for  mixed-sediment  layer  solids  as 
described  by  Chapra  and  Reckhow  (1983).  The  mass  balance  is 

0  -  v/A  -  (v,  *  v^4.(l  -  «p,  (8) 


The  framework  assumes  that  suspended  solids  concentration  sw  is  given. 
Therefore,  if  two  of  the  three  velocities  are  specified,  Equation  8  can  be 
employed  to  calculate  the  third. 

In  addition  to  these  transfer  velocities,  the  other  facet  of  the  model  related 
to  the  solids  budget  is  the  sediment  porosity.  Different  porosities  are  allowed 
for  the  mixed  layer  and  the  deep  sediments.  Within  these  regions,  porosity  is 
constant.  This  means  that  compaction  is  not  included. 
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Default  Parameters 


Several  parameters  in  the  above  framework  are  provided  by  the  software  as 
default  values.  These  are  suggested  values  that  represent  first  estimates.  They 
can  be  easily  modified  by  the  user. 


Volatilization  rate 

The  volatilization  rate  is  computed  by 


'w 


(9) 


where  vv  =  a  volatilization  transfer  coefficient,  m/year,  which  can  be  computed 
according  to  the  Whitman  two-film  theory  as 


v  = 


w 


K  H 
8  * 


(10) 


where  He~  a  dimensionless  Henry’s  coefficient  that  is  related  to  the  common 
form  He.  (atm  m3/gmole)  by 


H  = 

t 


(11) 


where 

R  =  universal  gas  constant  =  8.206  x  10'5  atm  m3/(gmole-kelvins) 

T  =  absolute  temperature,  kelvins. 

A  temperature  of  298K  (25  °C)  is  assumed  in  the  model. 

The  parameter  Kg  is  a  gas-film  mass  transfer  coefficient  (m/year)  that  can 
be  related  to  wind  speed,  Uw  (m/s),  and  molecular  weight,  MW  (g/gmole),  by 
(Mills  et  al.  1982) 


K  =  61320 
8 


18 


\).25 


MW 


U. 


(12) 


The  parameter  Kt  is  a  liquid-film  mass  transfer  coefficient  (m/year)  that  can 
be  related  to  wind  speed,  Uw  (m/s),  and  molecular  weight,  MW  (g/gmole),  by 
(Mills  et  al.  1982;  Banks  1975;  Banks  and  Herrera  1977) 
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(13) 


K, 


365 


32 


\>.25 


MW 


(0.728LC 


0.5 


-  0.317t/w  +  0.0372£/w2) 


Partition  coefficient 

The  partition  coefficient  for  organic  contaminants  is  computed  via 
(Karichhoff  et  al.  1979) 

K,  =  0.617/JC,.  (14) 

where 

foe  -  weight  fraction  of  organic  carbon  in  solid  matter,  g-orgC/g 
Km  =  octanol-water  partition  coefficient,  (mg/m3-octanol)/(mg/m3-water) 
Note  that  the  model  assumes  a  default  of  0.05  for  fM. 


Numerical  Methods 

The  coupled  set  of  differential  equations  is  solved  numerically.  The  two 
first-order  ODEs,  Equations  1  and  2,  are  solved  using  an  adaptive-step-size, 
fourth-order  Runge-Kutta  method  (Press  et  al.  1988).  The  partial  differential 
equation  is  solved  using  the  Crank-Nicolson  technique  (Chapra  and  Canale 
1988). 

The  total  number  of  years  for  which  the  model  is  run  is  determined  by 
approximating  the  time  required  for  the  contaminant  concentration  in  the  water 
to  decrease  to  10  percent  of  the  maximum  value  achieved.  This  is  determined 
by  using  an  analytical  solution  for  the  case  where  only  the  water  and  a  single 
well-mixed  layer  are  considered.  This  procedure  has  been  found  to  be 
effective  in  approximating  the  total  time. 
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3  Model  Confirmation 


The  time  variable  model  RECOVERY  has  been  applied  to  analyze  a  field- 
scale  experiment  in  which  a  flooded  limestone  quarry  was  dosed  with  equal 
quantities  of  the  insecticides  DDE  (dichloro-diphenyldichloro-ethylene)  and 
lindane  (hexachlorocyclohexane).  The  same  experiment  was  previously  ana¬ 
lyzed  with  a  time  variable  model  by  Di  Toro  and  Paquin  (1983)  with  good 
results. 


Study  Area 

The  study  was  conducted  in  two  flooded  limestone  quarries  located  near  the 
town  of  Oolitic  in  Bedford  County,  Indiana.  The  quarries  were  allowed  to 
flood  naturally  for  5  years  before  being  dosed  with  the  insecticides  lindane  and 
DDE.  One  quarry  was  experimentally  dosed  (quarry  T),  and  the  second  was 
used  as  the  control  (quarry  C)  (Waybrant  1973). 

Quarry  T  was  300  ft  long,  135  ft  wide,  and  an  average  of  50  ft  deep. 

Quarry  T  was  the  smaller  of  the  two  in  area,  with  extremely  clear  water  and 
generally  flat  bottom.  The  average  Secchi  disc  reading  for  quarries  C  and  T 
was  20  ft.  Both  quarries  exhibited  thermal  stratification  from  March  to 
November. 

The  relatively  soft  and  slightly  alkaline  water  found  in  the  quarries  tend  to 
indicate  that  the  quarries  received  very  little  groundwater  (Waybrant  1973). 

The  water  chemistry  and  physical  characteristics  of  the  quarries  indicated  that 
the  quarries  were  primarily  filled  with  rainwater  and  runoff  water  that  did  not 
percolate  through  the  limestone  to  the  groundwater. 

The  bottom  material  in  quarry  T  was  made  up  of  3  to  5  cm  of  fine  brown¬ 
ish  gyttja  underlaid  by  a  white  inorganic  mixture  of  limestone  dust  and  silica 
sand  (Waybrant  1973).  In  general,  the  top  sediment  layers  were  aerobic.  The 
mud  layer  contained  on  the  average  1  percent  sand,  42  percent  silt,  and 
57  percent  clay. 


Overview 


Quarry  T  was  treated  with  the  insecticides  lindane  and  DDE  at  a  concen¬ 
tration  of  200  parts  per  trillion  (ppt)  to  the  epilimnion  or  50  ppt  overall 
(2.7764  g  of  DDE  and  2.7752  g  of  lindane)  (Waybrant  1973).  The  quarries 
were  analyzed  after  the  treatment,  and  the  results  showed  that  essentially  all  of 
the  insecticides  were  initially  released  in  the  epilimnion.  The  quarries  were 
periodically  sampled,  and  the  results  are  presented  in  Table  2.  Both  water  and 
bottom  sediments  were  analyzed  for  the  insecticides  DDE  and  lindane.  The 
results  for  insecticides  in  the  bottom  sediment  material  are  shown  in  Table  3. 


Table  2 

Concentrations  of  DDE  and  Lindane  In  Water  Column 

Sampling  Day 

DDE,  ppt,  X  ±  SD 

Lindane,  ppt,  X  1  SD 

1 

44.0420  ±  16.336 

47.270 1  3.676 

5 

14.99413.429 

62.7921  10.419 

10 

4.270  ±  2.066 

60.828  1  2.733 

21 

3.5201  1.768 

58.8741  15.683 

42 

3.311  ±  1.324 

50.41812.087 

60 

2.13210.613 

29.8521  16.510 

81 

2.006  1  0.163 

32.939  1  4.048 

102 

1.46610.163 

34.577  1  4.059 

123 

0.988  1  0.220 

25.440  1  7.369 

144 

0.78910.171 

23.18912.378 

173 

0.958  10.163 

21.31810.667 

242 

0.95210.103 

20.284  1  1 .458 

Note:  SD  =  standard  deviation;  X  =  mean;  ppt  =  parts  per  trillion  (Waybrant  1973). 

The  concentration  of  DDE  and  lindane  remaining  in  the  water  column  after 
5  days  differed  significantly  (Table  2),  showing  different  transport  mechanisms 
by  the  two  insecticides.  After  Day  1,  a  significant  runoff  event  occurred  in 
quarry  T,  washing  a  significant  amount  of  sediment  into  the  quarry.  The  sus¬ 
pended  sediment  load  caused  the  DDE  concentration  in  the  water  column  to 
decrease  as  suspended  solids,  and  the  adsorbed  DDE  settled  to  the  lake  bottom. 
DDE  hydrophobic  properties  cause  affinity  to  suspended  solids  and  organic 
matter,  while  lindane,  a  more  polar  compound,  tends  to  have  higher  affinity  for 
water. 

Over  the  sampling  period,  an  extremely  rapid  distribution  of  DDE  in  the 
water  column  was  noted  as  opposed  to  an  87-percent  retention  of  lindane  in 
the  surface  layers  (Waybrant  1973).  The  quarry  underwent  turnover  at 
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Table  3 

DDE  and  Lindane  Concentration  in  Bottom  Sediments 

DDE,  ppb 

Lindane,  ppb 

Day 

Depth,  cm 

X  ±  SD 

%  of  Total 

X  ±  SD 

%  of  Total 

1 

1.5 

0.10 

7.05 

3.5 

7.05 

5.5 

5 

1.5 

18.66  ±  14.50 

30.27 

1.54  ±0.14 

2.50 

3.5 

5.5 

10 

1.5 

22.28  ±  6.04 

36.15 

1.48 

2.40 

3.5 

1.15±  1.49 

1.87 

<0.05 

5.5 

<0.10 

<0.05 

21 

1.5 

12.73  ±  11.30 

20.65 

1.64  ±0.43 

2.65 

3.5 

0.28  ±0.11 

0.45 

1.04  ±0.87 

5.5 

42 

1.5 

20.84  ±  5.53 

33.82 

0.98  ±  1.36 

1.59 

3.5 

0.14 

0.23 

5.5 

<0.10 

<0.05 

81 

1.5 

35.30  ±  27.45 

57.27 

1.13  ±  1.52 

1.83 

3.5 

<0.10 

0.99 

5.5 

<0.10 

0.95  ±  0.57 

123 

1.5 

27.27  ±  15.85 

45.24 

1.87  ±  1.77 

3.03 

3.5 

0.11 

0.18 

2.46  ±2.31 

3.99 

5.5 

0.88  ±1.06 

1.43 

1.40  ±0.97 

2.63 

173 

1.5 

30.92  ±  13.39 

50.17 

1.78  ±  1.60 

2.88 

3.5 

5.08 

8.24 

2.15  ± 

3.49 

5.5 

0.18 

<0.05 

242 

1.5 

19.69  ±2.80 

31.95 

<0.05 

3.5 

0.31  ±0.02 

0.50 

<0.05 

5.5 

0.39  ±0.45 

Note:  SD  =  standard  deviation;  X  =  mean;  ppb  =  parts  per  billion  (Way brant  1973). 

Day  144,  when  the  DDE  and  lindane  distribution  in  the  water  column  was 
essentially  homogeneous.  The  lindane  concentration  in  the  water  column 
decreases  until  Day  123,  after  which  the  insecticide  concentration  remains  con¬ 
stant  until  the  sample  after  the  following  spring. 

DDE  and  lindane  were  both  detected  in  the  bottom  sediments  at  signifi¬ 
cantly  higher  concentrations  after  the  initial  pesticide  loading.  However,  the 
ultimate  fate  of  both  pesticides  was  considerably  different.  DDE,  a  more 
hydrophobic  compound,  was  rapidly  deposited  in  the  bottom  sediments 
(Table  3),  while  lindane  concentration  never  reached  the  relatively  high  levels 
observed  with  DDE. 

DDE  peaked  at  35  parts  per  billion  (ppb)  in  the  bottom  sediments,  while 
lindane  levels  did  not  exceed  2.5  ppb  in  the  sediments.  Both  compounds 


reached  the  highest  level  in  the  sediments  between  Days  81  and  123,  lindane 
penetrating  the  lower  depth  of  the  bottom  sediments  while  DDE  remained 
mostly  on  the  top  layer. 

Some  important  notes  from  the  flooded  quarry  experiment  were  the  signif¬ 
icantly  higher  water  column  concentration  of  lindane,  an  order  of  magnitude 
above  DDE  at  the  end  of  the  study.  The  reverse  was  true  for  the  sediment 
concentration,  with  lindane  penetrating  to  the  deeper  sediment  layer  of  3.5  to 
5.5  cm.  In  addition,  lindane  profile  in  the  sediment  did  not  exhibit  a 
pronounced  vertical  gradient. 

Approximately  5  years  after  the  initial  dosing  of  the  quarry  on  June  21, 
1977,  several  sediment  samples  were  collected  and  analyzed  for  DDE  (Di  Toro 
and  Paquin  1983).  The  DDE  concentration  measured  in  that  sampling  event 
ranged  from  3.4  to  11.2  ppb  for  one  sample  and  from  2.9  to  4.2  ppb  for 
another.  The  water  column  was  sampled  at  the  same  time,  but  the  detection 
limit  of  the  analysis  was  30  ppt  (Di  Toro  and  Paquin  1983).  Sampling  was 
not  done  for  lindane  in  the  1977  sampling. 


Rate  Coefficients 

To  simulate  the  flooded  quarry  with  the  RECOVERY  model,  the  basic  as¬ 
sumptions  of  the  model  need  to  be  satisfied.  The  first  one,  the  water  body  is 
well  mixed,  might  be  satisfied  after  Day  23  (Table  2)  and  for  sure  at  Day  144. 
For  the  application  of  RECOVERY,  the  simulation  will  begin  between 
Days  81  and  123,  a  time  where  the  pesticide  profile  looks  fairly  uniform.  The 
second  one,  the  sediments  are  the  only  source  of  the  contaminant  to  the  water 
body,  is  satisfied  by  Day  81  for  DDE  and  Day  123  for  lindane.  At  those  days, 
the  concentration  in  the  sediment  peaks  for  each  of  the  simulated  compounds. 
Once  the  concentration  peaks  in  the  sediment,  the  water  cannot  be  a  source  of 
contaminant  anymore  (quasi-equilibrium  is  temporarily  established)  unless  a 
considerable  loss  of  contaminant  is  occurring  in  the  sediments. 

The  third  assumption,  the  initial  concentration  of  the  compound  in  the 
contaminated  region  is  uniform  throughout  that  region,  is  also  accomplished 
after  Days  81  and  123.  The  sediments  are  uniform  within  the  analytical  mea¬ 
surements,  and  in  the  deeper  sediments,  the  concentration  varies  in  the  vertical 
direction  only.  Therefore,  the  RECOVERY  simulations  will  start,  time  =  0,  at 
Day  81  for  DDE  and  Day  123  for  lindane. 

Initial  water  concentration  for  lindane  was  set  to  25.4  and  3.5  ppt  for  DDE. 
The  initial  sediment  concentration  was  set  to  35.3  ppb  for  DDE  and  1.87  ppb 
for  lindane.  Within  the  analytical  precision  of  the  data  and  the  collection 
scheme,  the  above  values  seem  to  be  representative  of  the  initial  conditions  for 
the  quarry  exposure.  The  inflow  to  the  quarry  was  estimated  to  be  insignifi¬ 
cant  in  comparison  to  the  overall  volume  of  water  in  the  quarry  (Di  Toro  and 
Paquin  1983).  The  water  depth  in  the  quarry  was  13.9  m,  and  the  quarry 
contained  5.23  x  107  {  of  water  (Waybrant  1973). 
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The  depth  of  the  mixed  sediment  layer  was  set  to  1.0  cm  for  DDE  and  3.0 
for  lindane  since  DDE  was  not  detected  at  substantial  levels  below  1.5  cm.  In 
this  case,  the  depth  of  the  deep  contaminated  sediments  was  set  to  the  allow¬ 
able  minimum,  which  was  0.05  m  (0.06  m  for  the  lindane  run)  below  which 
RECOVERY  can  become  unstable.  Suspended  solids  were  set  to  5  mg/C  based 
on  an  analysis  by  Di  Toro  and  Paquin  (1983)  on  Secchi  disc  readings  and  sedi¬ 
ment  trap  data.  Resuspension  was  set  to  zero,  the  burial  velocity  was  set  to 
5.0  x  104  m/year,  and  the  settling  velocity  was  estimated  from  the  solids  mass 
balance  in  RECOVERY.  The  estimated  value  was  87.5  m/year;  both  the  set¬ 
tling  velocity  and  the  burial  velocity  are  within  the  range  reported  by  Chapra 
(1983). 

The  partition  coefficient  for  DDE  was  computed  as  154,000  C/kg  based  on 
an  assumed  organic  content  of  0.05  (weight  fraction)  and  an  octanol/water 
partition  coefficient  of  5  x  104  (mg/m3  octanol)/(mg/m3  water).  The  molecular 
diffusivity  was  set  to  the  default  of  5  x  10'6  cm2/sec,  the  calculated  volatiliza¬ 
tion  rate  for  DDE  was  5.09  per  year,  and  the  degradation  rate  was  set  to  zero 
for  DDE  (half-life  of  15  years  (Howard  et  al.  1991)). 

For  lindane,  the  volatilization  rate  was  6.62  x  10'2  per  year,  and  the  degra¬ 
dation  rate  was  set  to  0.9  per  year  in  both  the  water  column  and  the  sediments. 
The  degradation  rate  agrees  with  the  value  of  0.0025/day  used  by  Di  Toro  and 
Paqin  (1983)  in  their  application.  A  detailed  listing  of  the  input  data  for  both 
the  lindane  and  DDE  runs  is  included  in  Appendix  B. 


Simulation  Results 


The  physical  and  chemical  parameters  estimated  and  selected  in  the  pre¬ 
vious  section  were  used  in  the  simulation  of  DDE  and  lindane  in  a  flooded 
limestone  quarry.  Results  from  the  initial  experiment  and  follow-up  work 
were  used  to  compare  with  the  simulation  from  RECOVERY. 

The  simulated  and  observed  lindane  distributions  are  shown  in  Figures  3 
and  4  for  water  and  bottom  sediments.  Lindane  decreases  from  25  ppt  to  less 
than  1  ppt  after  5  years  in  the  water  column.  Lindane  in  the  sediment  decreas¬ 
es  from  1.87  ppb  to  below  1  ppb  after  5  years.  Figure  5  shows  the  volatiliza¬ 
tion  of  lindane  from  the  water  column,  a  mechanism  in  the  depuration  of 
lindane  from  the  water  column.  Lindane’s  degradation  (0.9/year)  is  the  major 
mechanism  for  the  decrease  from  the  water  and  thus  the  flux  from  the  sedi¬ 
ments  back  into  the  water  column  (Figure  6). 

Figures  7  and  8  show  the  DDE  concentration  in  the  water  column  and  sedi¬ 
ment  bed,  respectively.  A  major  removal  mechanism  for  DDE  is  volatilization 
from  the  water  column  as  shown  in  Figure  9.  Initially,  DDE  was  transferred 
to  the  sediment  bed  by  deposition  of  suspended  sediment  contaminated  with 
DDE  from  the  quarry  initial  application.  After  the  initial  loading  to  the  bed, 
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DATA 


Figure  3.  Lindane  in  water  column 


Figure  4.  Lindane  in  sediment 


the  process  reverses  once  the  water  column  concentration  starts  decreasing 
because  of  the  loss  of  DDE  by  volatilization.  At  this  point,  the  sediment  bed 
becomes  the  source  of  DDE  for  the  water  column  (time  =  0,  see  Figure  10) 
since  there  is  no  degradation  in  the  bed;  the  only  removal  mechanism  is  vola¬ 
tilization.  Contrary  to  lindane,  where  the  pesticide  was  mostly  in  the  water 
column,  DDE  is  mostly  in  the  bed;  thus,  it  takes  longer  to  clean  the  bed  than  it 
does  to  clean  the  water.  If  degradation  in  the  sediment  bed  was  significant, 
then  a  much  faster  cleanup  would  occur. 
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Figure  5.  Volatilization  rate  of  lindane  from  water  computed  by  model 


{- 


Figure  6.  Flux  of  lindane  from  sediments  into  water  computed  by  model 

DDE  in  the  water  column  after  5  years  agrees  very  well  with  Di  Toro’s  (Di 
Toro  and  Paquin  1983)  prediction  of  less  than  1  ppt  (Figure  7).  In  the  sedi¬ 
ment,  the  simulation  (Figure  8)  shows  agreement  with  the  field  data,  but  a 
little  higher  than  the  Di  Toro’s  result  of  5.6  ppb  (Di  Toro  and  Paquin  1983). 
However,  both  models  predict  a  sediment  concentration  after  10  years  of 
approximately  1  ppb.  Figures  11  and  12  present  a  second  simulation  of  DDE 
in  the  water  and  sediment  bed,  respectively.  Porosity  in  the  mixed  bed  was 
increased  to  0.8,  and  the  degradation  of  DDE  in  the  water  was  set  to  3.0  per 
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Figure  8.  DDE  in  sediment 

year  to  account  for  photodegradation  of  DDE.  The  value  of  3.0/year  is  similar 
to  the  one  used  by  Di  Toro,  0.013/day  (4.75/year)  (Di  Toro  and  Paquin  1983). 
This  simulation  agrees  more  closely  with  the  other  model  results  and  shows 
remarkable  agreement  with  the  field  data  with  the  exception  of  spring 
sampling. 
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Figure  9.  Volatilization  rate  of  DDE  from  water  computed  by  model 


Figure  10.  Flux  of  DDE  from  sediments  into  water  computed  by  model 


Overall,  the  simulation  of  DDE  and  lindane  in  the  flooded  limestone  quarry 
can  be  simulated  with  confidence  using  the  model  RECOVERY.  The  confir¬ 
mation  application  was  performed  with  minimum  calibration.  In  addition,  very 
little  data  was  needed  to  accomplish  the  task.  There  are  three  very  important 
advantages  of  the  RECOVERY  model:  the  ease  of  application,  the  flexibility 
of  the  input  data  requirements,  and  the  user  friendly  environment  of  the  model. 
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Figure  1 1 .  Final  simulation  of  DDE  in  water  column 


Figure  12.  Final  simulation  of  DDE  in  sediment 
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4  User’s  Guide 


Use  of  the  RECOVERY  program  is  described  in  this  section.  Two  files  are 
required  to  run  this  program — RECOVERY.EXE  and  CORPCOMP.DAT. 
CORPCOMP.DAT  contains  the  chemical  database  information.  The  other  files 
included  on  the  distribution  disk  are  only  required  for  specific  hardware  sys¬ 
tems  and  to  alter  the  database.  (See  Changing  Chemical  Database  and  System 
Specific  Issues  below.)  Source  code  for  the  model  is  also  included.  (See 
Model  Source).  Do  not  run  RECOVERY  from  a  write-protected  disk,  as  this 
will  result  in  an  error. 

To  start  the  program,  move  to  the  directory  that  contains  the  software  and 
type  RECOVERY  at  the  DOS  prompt.  The  first  screen  is  a  title  screen.  Press 
ENTER  to  view  the  introduction  screen,  which  follows.  Press  ENTER  again, 
and  choose  between  creating  a  new  file  and  retrieving  data  from  an  existing 
data  file.  If  the  second  selection  is  made,  any  files  ending  in  .REC  and  the 
subdirectories  located  in  the  current  directory  are  displayed.  The  .REC  files 
are  the  data  files  that  are  created  when  saving  a  run  made  by  RECOVERY. 
(This  will  be  discussed  later).  If  your  file  is  not  located  in  the  current 
directory  or  in  one  of  the  listed  subdirectories,  press  BACKSPACE  and  specify 
the  path  name  to  your  file.  Your  .REC  file  should  then  be  displayed.  After  a 
file  is  chosen,  you  can  either  view  and/or  change  the  input  file  or  just  run  the 
model  using  that  data  file.  If  “create  a  new  file”  is  chosen,  a  list  of  com¬ 
pounds  is  displayed.  The  properties  of  18  compounds  are  included  in  the 
database.  If  the  compound  under  consideration  is  not  listed,  choose  number  0. 
The  selection  is  made  by  moving  the  cursor  to  the  desired  selection  using  the 
arrow  keys  and  pressing  ENTER. 

If  selection  0  was  made,  a  screen  appears  to  enter  the  name  of  the  contami¬ 
nant.  As  the  user  types,  the  name  will  appear  in  the  highlighted  area.  Correc¬ 
tions  can  be  made  by  using  the  BACKSPACE  or  DELETE  keys.  After  the 
name  has  been  typed  in,  press  ENTER  to  enter  the  name.  If  changes  are  nec¬ 
essary  at  this  point,  pressing  F2  will  allow  modification  of  the  name.  If  the 
user  starts  to  type  without  pressing  F2,  the  first  entry  is  lost  and  the  new  name 
appears.  This  editing  process  is  similar  to  that  used  by  Lotus  1-2-3  (R)  and  is 
used  on  all  input  screens.  Press  F10  to  store  the  data  and  advance  to  the  next 
screen. 
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The  next  five  screens  are  used  to  enter  or  change  input  values.  After  typ¬ 
ing  in  a  cell,  the  user  must  press  ENTER  before  proceeding  to  the  next  cell. 
The  first  morphometry  and  hydrology  screen  should  be  displayed.  If  a  differ¬ 
ent  compound  is  desired  at  this  point,  pressing  ESCAPE  will  allow  this 
change.  The  user  may  change  or  leave  the  water  concentration  at  the  default 
value  of  zero.  Three  of  the  next  four  values  need  to  be  entered.  Press  F9  to 
calculate  the  fourth,  and  its  values  will  be  displayed.  When  all  the  data  is 
entered,  press  F10  to  store  the  values  and  advance  to  the  next  screen.  If 
invalid  values  are  entered,  a  message  will  appear  at  the  bottom  of  the  screen. 

The  next  screen  is  for  the  contaminated  sediments  layer.  If  nothing  is 
changed,  the  default  values  will  be  used.  Press  F10  to  continue.  System  prop¬ 
erties  are  listed  on  the  next  screen.  Default  values  are  listed  for  the  first  eight 
variables.  These  values  can  be  changed  as  described  above.  Since  the  three 
velocity  values  are  interrelated,  two  of  the  three  values  need  to  be  entered. 
Typical  values  are  listed  in  the  center  column.  After  entering  two  values, 
press  F9  to  calculate  the  third  value.  It  will  appear  near  the  bottom  of  the 
screen.  Press  F10  to  store  the  data  and  proceed. 

The  next  two  screens  list  contaminant  properties.  Default  values  are  listed 
in  the  far  right  column  and  will  be  specific  to  the  selected  compound  if  it  is  in 
the  database.  These  values  will  be  used  unless  the  user  changed  them  in  the 
center  column.  The  second  screen  also  contains  intermediate  values  calculated 
from  the  previously  entered  data  (e.g.  ♦  Fdw  and  KjJ.  Press  F10  to  store  and 
advance. 

If  the  user  wants  to  return  to  a  previous  screen  in  the  input  process 
described  above,  he/she  can  press  ESCAPE  until  the  desired  screen  appears. 

The  next  screen  serves  as  a  check  to  see  if  the  entered  data  are  correct;  if 
not,  the  user  can  return  to  the  input  screens.  Selecting  the  third  option  starts 
the  model.  If  graphics  capabilities  are  available  on  the  system  being  used,  the 
concentration  profile  in  the  deep  sediments  is  displayed.  The  total  number  of 
years  is  calculated  as  the  approximate  time  when  the  water  concentration 
decreases  to  10  percent  of  its  maximum  value.  Since  the  numerical  method 
uses  an  adaptive-step-size,  the  time  between  profile  displays  will  not  be  even. 
When  the  computation  is  finished,  a  message  appears  in  the  lower  right-hand 
comer  of  the  screen. 

After  a  key  is  pressed,  another  menu  appears.  If  one  of  the  first  five 
options  is  chosen  or  number  8,  this  menu  appears  again  so  another  selection 
can  be  made.  Options  1-4  are  graphics  plots.  If  the  system  being  used  does 
not  have  graphics  capabilities,  a  message  appears  informing  the  user  of  this, 
and  a  key  can  be  pressed  to  return  to  the  menu.  Results  on  systems  having 
graphics  capabilities  other  than  those  listed  in  Chapter  4  of  this  report  will  be 
unpredictable.  The  first  two  plots  are  graphs  of  the  contaminant  concentration 
in  the  water  layer  and  the  mixed  layer,  respectively,  versus  time.  The  third 
plot  shows  the  flux  of  the  contaminant  from  the  sediments  to  the  water  versus 
time.  A  graph  of  the  volatilization  rate  of  the  contaminant  will  be  displayed  if 
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number  4  is  selected.  These  results  can  be  printed  by  typing  Shift+Print 
Screen.  Select  number  5  to  print  a  table  of  results  on  the  printer.  Number  6 
should  be  selected  to  make  modifications  to  the  existing  run,  and  number  7 
should  be  used  to  return  to  the  beginning  of  the  program.  Choose  number  8  to 
save  the  existing  run.  The  user  will  be  prompted  for  a  file  name  to  be  used  to 
store  the  results.  Choosing  number  9  will  result  in  exiting  the  program. 


Changing  Chemical  Database 

Currently,  data  for  18  chemicals  are  included  in  the  data  file, 
CORPCOMP.DAT A.  This  file  can  be  changed  and  enlarged  to  include  a 
maximum  of  39  compounds.  Changes  are  made  by  altering  the  file 
CORPSDAT.BAS,  which  is  included  with  the  distribution  disk.  This  file 
contains  the  following  data  for  a  compound: 

Name 

Molecular  Weight 
Octanol-Water  Partition  Coefficient 
Henry’s  Constant 
Molecular  Diffusivity 
Decay  Constants— water,  dissolved  phase 
water,  particulate  phase 
mixed  layer,  dissolved  phase 
mixed  layer,  particulate  phase 
deep  sediments,  dissolved  phase 
deep  sediments,  particulate  phase 

Additions  and  changes  are  made  by  altering  or  adding  data  statements  to  the 
file  CORPSDAT.BAS  and  running  the  program.  Errors  will  occur  in 
RECOVERY  if  there  are  more  than  39  compounds,  a  compound  name  with 
greater  than  28  letters,  or  additional  parameters  are  added  to  the  data 
statements. 


Model  Source 


In  order  to  access  the  model  source  code  in  the  QuickBasic  environment, 
type  qb  recovery  [\  reclib  when  you  are  in  a  directory  containing  QuickBASIC 
Version  4.5  and  the  following  files  (included  on  the  distribution  disk): 

RECOVERY.MAK 

RECOVERY.BAS 

CRANK.BAS 

ADAPTEST.BAS 

MODEL.BAS 

RECLIB.  QLB 


22 


Chapter  4  User’s  Guide 


Program  Features 


A  list  of  data  needed  to  run  this  program  can  be  found  in  Table  4.  Since 
the  compound  properties  for  several  contaminants  are  included  in  the  provided 
database,  the  user  needs  to  provide  compound  properties  only  if  the  contami¬ 
nant  under  consideration  is  not  one  that  is  in  the  database.  The  values  in  the 
database  that  comes  with  the  distribution  disk  are  listed  in  Table  5.  The  data 
file  can  be  changed  as  described  in  Chapter  5. 

The  data  are  entered  on  a  series  of  input  screens  with  editing  capabilities. 
The  user  can  easily  page  forward  and  backward  through  the  screens  while 
entering  data.  On-line  help  is  also  available.  If  any  of  the  entered  data  are 
invalid,  an  error  message  appears  instructing  the  user  to  re-enter  the  data.  For 
computers  with  graphic  display  capabilities,  the  deep-sediment  concentration 
profile  will  be  displayed  during  the  model  computation. 
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Table  4 

Input  Data  for  RECOVERY  Model 


Water  Column  Morphometry/Hydrology:  (3  specified,  4th  computed) 

Flow  through,  m3/year 
Residence  time,  year 
Surface  area,  m2 
Water  depth,  m 


Water  Layer: 

Initial  concentration,  jig/m3 


Contaminated  Sediments: 

Depth  of  contaminated  sediment  layer,  m 
Depth  of  mixed  layer,  m 
Mixed  sayer  surface  area,  m2 
Initial  concentrations,  |ig/m3 
Mixed  layer 

Deep  contaminated  sediments 


System  Properties: 

Suspended  solids  concentration  in  water,  g/m3 

Mixed  layer  porosity 

Deep  sediment  porosity 

Sediment  particle  density,  g/m3 

Wind  speed,  m/s 

Fraction  organic  carbon,  g-orgC/g-dry  wt  solids 
Water  layer 
Mixed  layer 

Deep  contaminated  sediments 

Two  of  the  following  three  velocities: 
Resuspension  velocity,  m/s 
Burial  velocity,  m/s 
Settling  velocity,  m/s 


Compound  Properties: 

Molecular  diffusivity,  cm2/s 
Decay  coefficients,  per  year:  particulate  -  water 
dissolved  -  water 
particulate  -  mixed  layer 
dissolved  -  mixed  layer 
particulate  -  deep  sediments 
dissolved  -  deep  sediments 
Henry’s  constant,  atm-m3/gmole 
Octanol-water  partition  coefficient 
Molecular  weight,  g/gmole 


Table  5 

Compound  Properties 

Compound 

Henry's  Constant 
@  25  °C 
atm  m3/gmole 

log  (Kj 

Molecular  Weight 
g/gmole 

Chlordane 

4.8  x  10* 

n 

2.78 

409.6 

DDT 

3.9  x  10* 

5.00 

354.5 

Dieldrin 

2.0  x  10* 

3.54 

381.0 

Lindane 

4.9  x  10* 

3.70 

290.0 

Arochlor  1 242 

8.4  x  10* 

6.72 

328.4 

Arochlor  1248 

3.5  x  10  3 

6.00 

295.0 

Arochlor  1254 

2.8  x  10* 

6.00 

326.0 

Arochlor  1 260 

7.1  x  10* 

6.50 

361.0 

Benzene 

5.5  x  10-3 

2.10 

78.0 

Chlorobenzene 

3.7  x  10  3 

2.98 

112.6 

Ethylbenzene 

8.7  x  10* 

3.15 

106.0 

Pentachlorophenol 

3.4  x  10* 

5.00 

266.0 

Phenol 

1.3  x  10* 

1.46 

94.1 

Toluene 

6.7  x  10* 

2.70 

92.0 

Anthracene 

1.7x10* 

4.45 

!  178.2 

Benzo(a)pyrene 

7.2  x  10* 

5.00 

252.0 

Naphthalene 

4.3  x  10* 

3.36 

128.0 

Chloroform 

4.2x10* 

1.95 

119.0 

Note:  Values  for  Henry's  constant  Kow  and  molecular  weight  are  from  Lyman  et  at.  (1982), 

Schnoor  et  al.  (1987),  and  Thomann  and  Mueller  (1987).  The  default  values  for  all 
molecular  diffusivities  are  5  x  10‘6  cm2/s.  The  default  values  for  all  decay  coefficients 
are  zero. 

The  final  results  are  displayed  as  plots  on  the  screen  and  can  be  printed  on 
a  printer  using  the  Shift+Print  Screen  command.  Numerical  results  can  be 
dumped  to  the  printer  also.  The  plots  are  as  follows: 

a.  Contaminant  concentration  in  the  water  versus  time. 

b.  Contaminant  concentration  in  the  mixed  layer  versus  time. 

c.  Flux  of  the  contaminant  from  the  sediments  into  the  water  versus  time. 

d.  Volatilization  losses  of  the  contaminant  from  the  water  versus  time. 
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The  input  data  can  be  stored  on  a  disk  and  retrieved  for  later  use.  The  user 
can  also  easily  return  to  the  input  screens  to  make  modifications.  Appendix  C 
gives  a  sequential  example  of  the  screens  in  running  the  RECOVERY  model. 

This  program  is  written  in  Microsoft  QuickBASIC  4.5  and  has  been  written 
for  the  following  combinations  of  graphic  adapters  and  displays: 

CGA  and  Color  Display 

EGA  and  Monochrome  Display 

EGA  and  Color  Display 

EGA  and  Enhanced  Color  Display 

VGA  and  Monochrome  Display 

VGA  and  Color  Display 

Hercules  Video  Card  and  Monochrome  Display 

A  math  coprocessor  is  not  required  to  run  RECOVERY.  The  use  of  one, 
however,  will  decrease  the  computation  time. 


System  Specific  Issues 

The  following  list  describes  some  procedures  that  are  specific  to  particular 
hardware  systems.  If  you  have  problems  running  RECOVERY,  please  consult 
the  list  below. 

a.  Use  of  Hercules  Graphics  Card 

If  your  system  has  a  Hercules  Graphics  Card,  you  will  need  to  load  the 
Hercules  driver  before  running  the  program.  The  driver, 
MSHERC.COM,  ((C)  Copyright  Microsoft  Corporation,  1987,  1988)  is 
provided  on  the  distribution  disk.  To  run  the  program,  type  MSHERC. 
If  you  do  not  run  this  program,  the  model  should  run  but  not  display 
any  graphics. 

b.  Strange  Sediment  Profile  Results 

If  your  system  has  all  of  the  following  characteristics,  you  may  have 
floating-point  math  problems. 

•  Uses  MS-DOS  Version  3.20. 

•  Boots  from  a  hard  drive. 

•  Has  a  math  coprocessor  (for  instance,  an  8087  chip). 

These  problems  can  be  avoided  by  booting  the  computer  from  a  floppy 

disk  drive  instead  of  the  hard  drive. 
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WARNING:  If  you  have  floating-point  math  problems,  the  model 
information  will  not  be  correct.  This  should  be  evident  by  the 
graphical  results. 

c.  Use  of  a  Compaq  386  with  a  Math  Coprocessor 

If  you  are  using  a  Compaq  386  that  has  a  math  coprocessor  installed, 
RECOVERY  may  crash.  In  order  to  use  the  model,  turn  the  coproces¬ 
sor  off  by  adding  the  line  SET  N087  =  ’  ’  to  the  CONFIG.SYS  file 
and  rebooting  your  machine. 
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5  Conclusions  and 
Recommendations 


RECOVERY  is  designed  for  interactive  implementation  via  a  personal  com¬ 
puter.  The  program  allows  the  user  to  rapidly  generate  and  analyze  recovery 
scenarios  for  contaminated  sediments.  The  software  includes  graphical  dis¬ 
plays  and  is  self-documented.  The  intent  of  this  report  is  to  document  the 
model  and  confirm  it  against  existing  field  data. 

The  confirmation  of  RECOVERY  against  the  quarry  experiments  shows  the 
ability  of  the  program  to  simulate  behavior  of  organics  in  a  real  system  with  a 
limited  amount  of  data.  Results  are  similar  to  those  predicted  by  another 
model  and,  thus,  appear  to  demonstrate  the  validity  of  the  algorithms  used  to 
describe  the  fate  and  transport  of  organics. 

One  problem  found  in  the  confirmation  exercise  was  the  simulation  of  the 
initial  part  of  the  quarry  experiment  where  an  initial  organic  load  was  applied 
to  the  pond.  The  current  formulation  of  the  model  assumes  no  loadings  asso¬ 
ciated  with  the  inflow.  To  overcome  the  loading  problem  and  the  completely 
mixed  assumption,  the  simulation  was  started  at  a  time  where  the  sediments 
were  contaminated  with  pesticides,  and  the  surface  water  pesticide  concentra¬ 
tion  was  fairly  uniform  throughout  the  pond.  After  satisfying  the  model  as¬ 
sumptions,  the  simulation  proceeded  with  good  results  when  compared  with 
both  the  data  collected  at  the  site  and  another  model. 

The  RECOVERY  model  can  be  applied  in  a  variety  of  scenarios,  including 
lakes,  embayments,  harbors,  estuaries,  and  ocean  parcels,  as  long  as  the  as¬ 
sumptions  of  a  completely  mixed  water  body  is  acceptable  along  with  other 
assumptions  listed  in  Table  1.  Average  yearly  flushing  would  be  needed  in 
any  of  the  applications  above,  in  addition  to  the  other  chemical  and  physical 
parameters  for  a  well-defined  simulation. 

Mechanistic  and/or  modular  improvements  recommended  for  RECOVERY 
are  the  following: 

a.  Allow  external  contaminant  loadings  from  point  and  nonpoint  sources. 

b.  Allow  atmospheric  contaminant  loads,  such  as  dry  or  wet  deposition. 
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c.  Allow  the  variation  of  physical  and  chemical  properties  in  the  contami 
nated  sediment 

d.  Include  food  chain  uptake,  such  as  fish  gill  exchange  and  sediment 
uptake. 

e.  Develop  a  version  of  RECOVERY  to  assess  the  impact  on  surface 
waters  from  bottom  sediments  contaminated  with  heavy  metals. 

The  above  recommendations  include  both  short-term  modifications  (the  first 
two)  and  long-term  improvements.  The  short-term  modifications  will  provide 
a  wider  range  of  applications  for  the  RECOVERY  model  and  would  not  be 
major  modifications  of  the  program.  Recommendations  c  and  d  are  a  major 
undertaking  and  would  imply  substantial  modifications  of  the  program.  The 
last  recommendation  would  involve  substantial  effort  since  both  transport  and 
kinetic  algorithms  would  need  modifications  and/or  development. 
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Appendix  A 
Nomenclature 


Am  Mixed  layer  surface  area,  m2. 

Aw  Water  surface  area,  m2. 

C,  Contaminant  concentration  in  the  inlet  to  the  water  body, 

pg/m3. 

Cm  Contaminant  concentration  in  the  mixed  layer,  pg/m3. 

CmO  Initial  contaminant  concentration  in  the  mixed  layer,  pg/m3. 

Cs  Contaminant  concentration  in  the  deep  sediments,  pg/m3. 

Initial  contaminant  concentration  in  the  deep  sediments,  pg/m3. 

C/O)  Contaminant  concentration  at  the  top  of  the  deep  sediments, 

pg/m3. 

Cw  Contaminant  concentration  in  the  water,  pg/m3. 

Cw0  Initial  contaminant  concentration  in  the  water,  pg/m3. 

Ds  Molecular  diffusivity,  m2/s. 

Dm  Diffusion  coefficient  in  the  sediment  pore  water,  m2/s. 

Ratio  of  contaminant  concentration  in  the  mixed  layer  pore 
water  to  concentration  in  the  total  mixed  layer,  dimensionless. 

Ffa  Ratio  of  contaminant  concentration  in  the  deep  sediment  pore 

water  to  concentration  in  the  total  deep  sediments, 
dimensionless. 

F^  Fraction  of  contaminant  in  particulate  form  in  the  water, 

dimensionless. 
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foe 

H' 

H, 

J 

J ms 

K* 


K, 


Kn, 


ky 

kyy 

MW 

Q 

R 

Syy 

t 

T 

\ 


Weight  fraction  of  organic  carbon  in  the  solid  matter, 
g-orgC/g-solid. 

Fraction  of  contaminant  in  dissolved  form  in  the  water, 
dimensionless. 

Henry’s  law  constant,  dimensionless. 

Henry’s  law  constant,  atm-m3/gmole. 

Flux  of  the  contaminant,  g/(m2-year). 

Flux  of  the  contaminant  from  the  mixed  layer  to  the  deep 
sediments,  g/(m2-year). 

Contaminant  partition  coefficient  in  the  water,  m3/g. 

Contaminant  partition  coefficient  in  the  sediments,  m3/g. 

Gas-film  mass  transfer  coefficient,  m/year. 

Liquid-film  mass  transfer  coefficient,  m/year. 

Decay  rate  constant  of  the  contaminant  in  the  mixed  layer, 
/year. 

Octanol-water  partition  coefficient,  (mg/m3-Octanol)/ 
(mg/m3-water). 

Decay  rate  constant  of  the  contaminant  in  the  deep  sediments, 
/year. 

Volatilization  rate  of  the  contaminant,  /year. 

Decay  rate  constant  of  the  contaminant  in  the  water,  /year. 
Molecular  weight  of  the  contaminant,  g/gmole. 

Flow  rate  into  and  out  of  the  water  body,  m3. 

Universal  gas  constant,  atm-m3/(gmole-kelvins). 

Suspended  solids  concentration  in  the  water,  g/m3. 

Time,  years. 

Absolute  temperature,  kelvins. 

Burial  velocity,  m/year. 
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\d  Diffusive  mass  transfer  coefficient,  m/year. 

Vm  Volume  of  the  mixed  layer,  m3. 

vr  Resuspension  velocity  of  the  sediments  into  the  water,  m/year. 

v.  Settling  velocity  of  the  particulate  matter  from  the  water  to  the 

sediments,  m/year. 

vv  Volatilization  transfer  coefficient,  m/year. 

Vw  Volume  of  the  water,  m3. 

Uw  Wind  speed,  m/year. 

z’  Length  over  which  the  concentration  gradient  exists  at  the 

sediment/water  interface,  m. 

zm  Depth  of  the  mixed  layer,  m. 

zw  Depth  of  the  water  layer,  m. 

p  Density  of  the  sediment  solids,  g/m3. 

<|>  Porosity 
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Appendix  B  Input  Screens  for  DDE  and  Lindane  Runs 


RECOVERY  RUN 
DDE 


l 


Hater 


Hixed 


Deep 
|  Contaminated 
Sediments 


Clean 

Sediments 


Heip-Fl 


HGRPHOflETRY  AMD  HYDROLOGY 


COfiTfinWAIED  SEDMEMIS  LAYEB 


L  =  Contaminated  Sediment  Depth  (m) 

2  =  Depth  of  fiixed  layer  (cm) 
Hixed  Layer  Surface  Area  (n*2) 

Initial  Concentrations  (Hg/n*3) 
Hixed  Layer: 
Deep  Contaminated  Sediments: 


iB09E*ea 


3.5B0E+A4 

3.500E*A4 


Edit-F2 


Store  Data  and  Aduance-F16 


INPUT  SCREEN  2 


RECOVERY  RUN 
DDE 


SVSIB1  PROPERTIES 

Suspended  Solids  Concentration  in  Hater  (3/11*3) 
Kixed  Sedinent  Porosity 

B.65 

Deep  Sedinent  Porosity 

B.65 

nixed  Sedinent  Particle  Density  (g/n*3) 

2.58EH96 

Deep  Sedinent  Particle  Density  (g/n*3) 

2,5BE<86 

Hind  Speed  (Vsec) 

2.88 

Height  fraction  carbon  in  solid  -  water 

5.88E-82 

Height  fraction  carbon  in  solid  -  nixed  layer 

S.B0E-B2 

Height  fraction  carbon  in  solid  -  sedinents 

5.BBE-82 

INPUT  TUO  OF  THE  FOLLOWING  VELOCITIES  AND  PRESS  F9= 

Velocity  (n/yr)  1  Typical  Value 

I  Value  for  this  Run 

8,eeE<ee 

S.08E-84 


8.75E*B1 


INPUT  SCREEN  3 


DDT  PROPERTIES 

Property 

I  Value  for  This  Run  1 

Default  Value 

nry's  Con,  (atn-n j/gnole) 
lecular  Height 

Dctanol-Hater  Partition  Coeff . 
(ng/n*3  octanolVtng/n^  water) 


354.58 

5.00E*84 


9.92E-81 

7.6SE-83 

1.73E-83 

7.41E-84 


1.54E-03 

1.54E-03 

1.54E-83 

5.09E*00 


INPUT  SCREEN  4 
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3.98E-85 

354.58 

5.80EH94 


9.92E-01 


RECOVERY  RUN 2 
DDE 


NORPHONETRY  AND  HYDROLOGY 


Initial  Concentration  of 
DDI  0*5/11*3) : 


Deep 

Contaninated 

Sedinents 

Clean 

Sedinents 


BITER  3  OF  THE  FOLLOWING  4  UAIJUES  AMD  PRESS  F* 


Water  Surface  Area  (n*2) 
Water  Depth  (n) 
Flog  Through  (n*3/yr) 
Residence  Tine  (years) 


Calculated  Flow  Through  (n*3/yr): 


Edlt-F2  ReCalc-F9  Store  Data/Ad wnce-F  10 


Deep  L 

Contaninated 
Sedinents 

Clean 

Sedinents 


NORPHONETRY  AND  HYDROLOGY 


CONTANINATED  SEDINENTS  LAYER 


L  =  Contaninated  Sedinent  Depth  (n): 
2  =  Depth  of  Nixed  Layer  (cn): 
Nixed  Layer  Surface  Area  (n*2): 

Initial  Concentrations  0*g/n*3) 
Nixed  Layer: 
Deep  Contaninated  Sedinents : 


1.09 

LB89E+G0 


3.500E+R4 

3.500EHJ4 


Store  Data  and  Aduance-F10 
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RECOVERY  RUN 2 
DDE 


SYSIH1  PROPERTIES 

Suspended  Solids  Concentration  in  Uater  (g/*3) 

■m 

nixed  Sedinent  Porosity 

B.80 

Deep  Sedinent  Porosity 

0.05 

nixed  Sedinent  Particle  Density 

2,S8E*B6 

Deep  Sedinent  Particle  Density  (g/n*3) 

2.5BE«B6 

Hind  Speed  (n/sec) 

2,88 

Height  fraction  carhon  in  solid  •  uater 

S.BBE-B2 

Height  fraction  carbon  in  solid  -  nixed  layer 

5.88E-82 

Height  fraction  carbon  in  solid  -  sedlnents 

S.00E-B2 

INPUT  TWO  OF  THE  FOLLOWING  VELOCITIES  MD  PRESS  F9= 

Uelxity  Ot/yr)  |  Typical  Value 

1  Value  for  this  Run 

e.eeE>89 

S.0BE-04 


B.7SE«B1 


INPUT  SCREEN  3 


DDT  PROPERTIES 

Property 

1  Value  for  This  Run  1  Default  Value 

lecuiar  Diffusiuity  (cn  2/sec) 


ECAY  C0EFF1CIBITS  tl/yr) 

issolved  CcmtaninenV 

In  Water 
In  nixed  Layer 
In  Deep  Sedlnents 

articulate  Contaninenf- 

In  (later 
In  nixed  Layer 
In  Deep  Sedlnents 


Store  Data  and  Aduance-FIB 


INPUT  SCREEN  4 
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RECOVERY  RUN 
LINDANE 


ttORPHOfiETRV  m  HVDBOLOGV 


Initial  Concentration  of 

— >  Lindane  Oym*3): 

Hater 

Hixed 

BITER  3  OF  THE  FOLLOWING  4  VALUES  AW)  PRESS  F9= 

Jeep 

Hater  Surface  Area  Ort): 

1800 E +00 

Contaminated 

Hater  Depth  (m): 

10.8 

Sediments 

Flou  Ihrough  (m*3/yr): 

Residence  Time  (years): 

58.00 

Clean 

Sediments 

Calculated  Flow  Through  (n*3^r):  280.80 

-FI  Edlt-F2  ReCalc-F9  Store  Data/Aduance-F10  Change 

Conpound-Esc 

ttORPHOftElRV 


HVDROLOGV 


Deep  L 

Contaminated 
Sediments  | 

Clean 

Sediments 


COrtTArt  MATED  SEDMENTS  LAVER 


L  =  Contaminated  Sediment  Depth  (m): 
2  *  Depth  of  nixed  Layer  (cm): 
Hixed  Layer  Surface  Area  ("*2): 

Initial  Concentrations  (pg/n*3) 
nixed  Layer: 
Deep  Contaminated  Sediments: 


3.88 

1B08E+90 


1.879E+03 

2.408E+83 


Stare  Data  and  Aduance-FIB 
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RECOVERY  RUN 
LINDANE 


SYSIEJI  PROPERTIES 

Suspended  Solids  Concentration  in  Uater  (g/n*3) 

nixed  Sedinent  Porosity 

B.65 

Deep  Sedinent  Porosity 

8.6S 

nixed  Sedinent  Particle  Density  (g/n*3) 

2.5BE<86 

Deep  Sedinent  Particle  Density  (g/n*3) 

2.SBE«86 

Wind  Speed  (n/sec) 

1.04 

Weight  fraction  carbon  in  solid  -  uater 

S.B0E-B2 

Height  fraction  carbon  in  solid  -  nixed  layer 

5.BBE-62 

Height  fraction  carbon  in  solid  -  sedinents 

5.B0E-B2 

INPUT  TWO  OF  THE  FOUOUIHG  VELOCITIES  AND  PRESS  W 

Velocity  (Vyr)  1  Typical  Value 

I  Value  for  this  Run 

B.8BE»88 

5.34E-04 

1.0BE+Q2 


INPUT  SCREEN  3 


UNDONE  PROPERTIES 

Property 

1  Value  for  This  Run  1  Typical  Value 

blecular  Diffusivity  (cn*2/sec) 


5.0E-86 


ECAY  C0ETF1CIBITS  (1/yr) 

issolued  Contaninerrt : 

In  Uater 
In  Rixed  Layer 
In  Deep  Sedinents 

articulate  Contaninent: 

In  Uater 
in  Hired  Layer 
In  Deep  Sedinents 


Store  Data  and  Aduance-F18 


INPUT  SCREEN  4 
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RECOVERY  RUN 
LINDANE 


UHDf)NE  PROPERTIES 

Property 

Value  for  This  Run 

Default  Value 

fenry's  Con.  (atn-n^gnole) 

4.98E-07 

lolecular  Height 

298.80 

298.08 

tetanol-Vater  Partition  Coeff. 
(ngV3  octanolV(ng4f3  uater) 

5.B1E+03 

5.01E+83 

Fdu 

9.99E-81 

9.99E-01 

Fpu 

7.72E-84 

7.72E-04 

Fdpn 

7.36E-03 

7.36E-83 

"dps 

7.36E-83 

7.36E-03 

1.55E-84 

1.55E-04 

1.55E-04 

l.SSE-84 

1.55E-84 

1.55E-B4 

iji  pm  ipyiTp)ip7!t%yi^iy^tqv; 

■  ■■  B 

6.62E-82 

6.62E-82 

ReCalc-F9  Store  Data  and  Adgance-F10 

INPUT  SCREEN  5 


Appendix  B  Input  Screens  for  DDE  and  Lindane  Runs 
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Appendix  C  Screens 


ci 


RECOVERY 

A  Decision  Support  Model  for  Screening  the  Fate  of  In-Place 
Contaminated  Sediments  in  Aquatic  Environments 


June.  1989 

Beta  Test  Vers  ion  l . 0 

Written  for  The  U.S.  Army  Corps  of  Engineers 
Waterways  Experiment  Station 
and 

New  Eng i and  Division 

Written  by  Jean  M.  Boyer  and  Steven  C.  Chapra 

Center  for  Advanced  Decision  Support 
in  Water  and  Environmental  Systems 
University  of  Colorado 
Boulder,  CO  80309-0428 

MsHerc  program  (C)  Copyright  '87,  *88  Microsoft  Corp.  All  rights  reserved. 


Press  any  key  to  continue. 


SCREEN  1 


This  program  is  a  tool  for  screening  the  fate  of  tn-place  contaminated 
sediments  In  aquatic  :n v i r onmen ts .  It  predicts  the  concentration  of  a 
contaminant  In  the  water,  the  mixed  sediment  layer  and  in  the  deep  sedi¬ 
ments  over  time.  The  flux  of  the  contaminant  from  the  sediments  into  the 
water  is  also  predicted.  Note  that  the  total  number  of  years  for  which 
the  model  Is  run  is  determined  by  approximating  the  time  required  for  the 
toxic  concentration  in  the  water  to  decrease  to  10%  of  the  maximum  value 
achieved,  up  to  a  maximum  of  100  years. 


Please  refer  to  the  documentation  accompanying  this  software  for  de¬ 
tailed  instructions  and  model  theory. 


(f  you  need  help  while  entering  data  into  an  input  screen,  press  FI. 

Press  Ctrl*Break  to  terminate  the  program  at  any  time.  Press  ShifWrint 
Screen  to  print  a  graph  displayed  on  the  screen. 


Press  any  key  to  continue. 


SCREEN  2 
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Highlight  your  choice  and  press  Enter  or 
press  the  number  preceding  your  selection. 


1.  Create  a  new  file. 

2.  Retrieve  data  from  an  existing  file. 


These  ate  the  .REC  files  and  subdirectories  listed  in  the  specified 
directory.  Highlight  one  of  these  using  the  up  and  down  arrows  and 
press  Enter  or  press  Backspace  to  specify  an  alternative  directory. 
Press  escape  to  return  to  the  previous  screen. 


Directory:  C:\QB45\ 


CHLOR.REC 


ANTH.REC 
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The  properties  of  the  following  contaminants  are  included  in  the  database. 
If  the  contaminant  you  want  to  consider  is  not  listed,  choose  0. 


Highlight  your  choice  using  the  arrow  keys  and  press  ENTER. 


0. 

Compound  Not  Listed 

10. 

Chlorobenzene 

1. 

Chlordane 

11. 

Ethyl  Benzene 

2. 

DDT 

12. 

Pentachlorophenol 

3. 

Dieldrin 

13. 

Phenol 

4. 

Lindane 

14. 

Toluene 

5. 

Arochlor1242 

15. 

Anthracene 

6. 

Arochlor1248 

16. 

Benzo(a)pyrene 

7. 

Arochlor1254 

17. 

Naphthalene 

8. 

Arochlor1260 

18. 

Chloroform 

9. 

Benzene 

SCREEN  4b 


Highlight  your  choice  and  press  Enter  or 
press  the  number  preceding  your  selection. 


1 .  View/Change  Input  file. 


2.  Run  Using  Input  File. 


SCREEN  5 
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1 

CHLORDANE  PROPERTIES 

\  Property 

Value  for  This  Run 

Default  Value 

Henry's  Con.  (atm-mA3/gmole) 

4.80E-05 

4.80E-05 

Molecular  Weight 

409.60 

409.60 

Octanol-Water  Partition  Coeff. 

6.03E+02 

6.03E+02 

(mg/mA3  octanol)/(mg/mA3  water) 

CALCULATED  DATA 

Fdw 

1.00E+00 

1.00E+00 

| 

Fpw 

3.72E-05 

3.72E-05 

I 

Fdpm 

1.28E-01 

1J28E-01 

■ 

Fdps 

1.28E-01 

1.28E-01 

i 

Partition  Coefficients 

Kdw 

1.86E-05 

1.86E-05 

Kdm 

1.86E-0.5 

1.86E-0.5 

Kds 

1.86E-0.5 

1.86E-0.5 

■ 

Volatilization  rate  (1/yr) 

1.14E+01 

1.14E+01 

i 

Help-FI  Edit-F2 

1 

ReCaIc-F9  Store  Data  and  Advance-FI  0 

1 

Highlight  your  choice  and  press  Enter  or 
press  the  number  preceding  your  selection. 


1 .  Return  to  input  screens.  Change  compound. 


2.  Return  to  input  screens.  Same  compound. 


3.  Entered  data  are  ok.  Continue. 
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Concentration  Profile  Below  the  Mixed  Sediment  Layer  (depth  in  m) 
Chlordane 


D 

e  1.0 
P 


10  20  30  40  SO  60  70  80  90  100 


Concentration  (micrograms/mA3)  x  10A1 


Press  any  key  for  output 


Total  years:  1.0 


Highlight  your  choice  and  press  Enter  or 
press  the  number  preceding  your  selection. 


1 .  Plot  Contaminant  Concentration  in  Water. 

2.  Plot  Contaminant  Concentration  in  Mixed  Layer. 

3.  Plot  Flux  of  Contaminant  from  Sediments  into  Water. 

4.  Plot  Volatilization  Rate  from  the  Water. 

5.  Print  Table  of  Results  to  Printer. 

6.  Return  to  Input  Screens. 

7.  Return  to  Beginning  of  Program. 

8.  Save  the  Input  Data  for  This  Run. 

9.  Exit  Program. 
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Chlordane  Concentration  in  the  Water  (micrograms/mA3) 


C 

o 

n 

c 

e 

n 

t 

r 

a 

t 

i 

o 

n 


Time  (Years) 


Press  Shift+PrintScreen  to  print. 


Press  any  key  to  continue. 


SCREEN  14 
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GC  CO 
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